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ABSTRACT 



We present new multi-wavelength millimeter interferometric observations of the Herbig Ae star HD 163296 obtained with the 
IRAM/PBI, SMA and VLA arrays both in continuum and in the ^'^CO, ^^CO and C^^O emission lines. Gas and dust properties 
have been obtained comparing the observations with self-consistent disk models for the dust and CO emission. The circumstellar 
disk is resolved both in the continuum and in CO. We find strong evidence that the circumstellar material is in Keplerian rotation 
^ ' around a central star of 2.6 Mq. The disk inclination with respect to the line of sight is 46°±4° with a position angle of 128°±4°. 
I The slope of the dust opacity measured between 0.87 and 7 mm {f3 = 1) confirms the presence of mm/cm-size grains in the disk 
midplane. The dust continuum emission is asymmetric and confined inside a radius of 200 AU while the CO emission extends 
up to 540 AU. The comparison between dust and CO temperature indicates that CO is present only in the disk interior. Finally, 
, we obtain an increasing depletion of CO isotopomers from ^^CO to ^''CO and C^*0. We argue that these results support the 
■ idea that the disk of HD 163296 is strongly evolved. In particular, we suggest that there is a strong depletion of dust relative 
' to gas outside 200 AU; this may be due to the inward migration of large bodies that form in the outer disk or to clearing of a 
C . large gap in the dust distribution by a low mass companion. 
O 



1. Introduction 

Millimeter and sub-millimeter interferometers are provid- 
ing an increasingly detailed description of disks around 
pre-main sequence stars of solar (T Tauri stars; TTS) 
and intermediate mass (Herbig Ae; HAe). Both dust con- 
tinuum emission and emission in molecular lines are ob- 
served and spatially resolved in a number of disks, yield- 
ing information on the disk density and temperature, the 
dust properties and the gas chemistry and dynamics in the 
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outer disk (e.g., Natta et al [20071 Dutrey et al. [20061 and 
references therein). Thanks to the recent instrumental im- 
provements, it is now possible to build upon the original 
detections and study more accurately the disk structure 
details. The number of well-studied disks is however still 
very small, practically restricted to the most massive and 
luminous ones; still, it is clear that disks differ from one 
another. Recently, it has been reported evidence of spiral 
structures in AB Aur, a 2-3 Myr old intermediate mass 
star, and of deviations from Keplerian rotation (Pietu et 
al. mm Corder et al. 1^0051) : the classical TTS LkCal5 
has a large inner hole of size ~ 50 AU depleted of dust, 
while the HAe star MWC 480 has a smooth disk with 
an optically thick (at millimeter wavelengths) inner re- 
gion of radius ^ 35 AU (Pietu et al. I2006p . Both spiral 
structures and large gaps are evidence of dynamical per- 
turbations, possibly due to the effect of large planets. The 
existence of both unperturbed and distorted disks among 
pre-main sequence stars suggests that the planet forma- 
tion is actively occurring during this evolutionary stage, 
leaving detectable marks on the parent disks. It is there- 
fore important to study in detail as many disks as possible, 
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in order to characterize their basic properties and to de- 
tect deviations from the simple patterns of homogeneous 
disks in Keplerian rotation. 

We report in this paper a detailed study of the disk 
associated to the HAe star HD 163296, using observa- 
tions in the continuum and CO lines obtained with three 
different interferometers, namely the Very Large Array 
(VLA) at 7mm, IRAM Plateau de Bure Interferometer 
(PBI) at 1.3 and 2.6 mm and the Submillimeter Array 
(SMA) at 0.87 mm. HD 163296 is a star of spectral type 
Al, mass of roughly 2.3 Mq, distance 122 pc (van den 
Ancker et al. L19_98!). Early OVRO observations (Mannings 
and Sargent 1997, 1 have shown the presence of a disk with 
a minimum mass ~ 0.03 Mq and evidence of rotation from 
the CO lines. The disk is seen in scattered light by Grady 
et al. (j2000l [T999|) . with radius of ~ 500 AU; it has an 
associated jet seen in Ly-a with HST, extending on both 
sides of the disk orthogonally to it (Devine et al. 120001 
Wassell et al. [2006l) . Natta et al. ((2003)) found evidence of 
evolved dust in the outer disk of HD 163296 by comparing 
the VLA 7 mm flux to the OVRO observations. 

The results we present here have much higher spatial 
resolution and wavelength coverage than what has been 
previously reported. They allow us to measure accurately 
the dynamics of the disk as well as the disk and dust 
properties and to test the capability of disk models to 
account for the observations. As we will show, they suggest 
that the HD 163296 system is probably evolving towards 
a debris disk phase. 

The structure of the paper is as follows. Sec. [H will 
describe the observations. The results will be presented in 
Sec.[3l where we will derive some of the disk parameters. A 
more detailed analysis, using self-consistent disk models of 
the dust and CO line emission will be presented in Sec. [U 
Sec. [5] contains the results, which will be further discussed 
in Sec. [S] Summary and conclusions follow in Sec. [7] and 
1 



2. Observations and data reduction 

2.1. PBI observations 

The PBI observations were carried over the 2003/2004 
winter season. The six 15 m dishes were used in the most 
extended configuration providing a baseline coverage be- 
tween 25 and 400 m. The corresponding angular resolu- 
tions are reported in Tab. [TJ The receivers were tuned to 
observe the ^^CO J=2-l line and the nearby continuum at 
1.3 mm, while at 2.8 mm the ^^CO J=l-0, and C^^O J=l- 
lines were observed along with the continuum. Bandpass 
and complex gain calibrations were ensured by observa- 
tions of standard IRAM calibrators. The phase stability 
was excellent during our observations and only a minimal 
amount of editing of the data was necessary. All calibra- 
tions were performed using the standard CLIC suite of 
programmes within the GILDAS software package. The 
calibrated uv data were then exported for the subsequent 



analysis. The accuracy of the flux density scale calibration 
is expected to be within 20% at these wavelengths. 

2.2. SMA observations 

The SMA observations of HD 163296 were made on 
August 23rd, 2005 using the Compact Conflguration of 
seven of the 6 meter diameter antennas, which provided 
21 independent baselines ranging in length from 8 to 80 
meters. The SMA digital correlator was conflgured with 
a narrow band of 512 channels over 104 MHz, which pro- 
vided 0.2 MHz frequency resolution, or 0.18 km s~^ veloc- 
ity resolution at 345 GHz, and the full correlator band- 
width was 2 GHz. The weather was good with r(225 
GHz) around 0.06 and the double-sideband (DSB) sys- 
tem temperature were between 200 and 500 K. The source 
HD 163296 was observed from HA -3 to 4.5. Calibration of 
the visibility phases and amplitudes was achieved with ob- 
servations of the quasar 1921-293, typically at intervals of 
25 minutes. Observations of Uranus provided the absolute 
scale for the flux density calibration and the uncertainties 
in the flux scale are estimated to be 20%. The data were 
calibrated using the MIR software packag^ll- 

2.3. VLA observations 

HD 163296 was observed at the NRAO/VLA as part of 
a larger survey for 7 mm disk emission around Herbig Ae 
stars (see Natta et al. . 120041) . Data were obtained with the 
array in the C and D configurations, in several occasions 
from Dec 2001 through May 2003. Accurate pointing was 
ensured by hourly pointing sessions at 3.6 cm on a bright 
extra galactic object. The array offered baselines from the 
shadowing limit through ~ 3.4 km, although some data 
was obtained in a hybrid DnA configuration, all the data 
from the longer baselines had to be rejected due to large 
phase fluctuations that could not be corrected. The result- 
ing uv-plane coverage offered an angular resolution of 0.5" 
and a maximum recoverable size of the order of ~40". All 
the data was edited and calibrated using standard recipes 
in AIPS. The short term complex gain variations were 
corrected using frequent (few minutes cycle) observations 
of the quasar 1820—254, while the flux scale was set ob- 
serving the VLA calibrator 1331-1-305. This procedure is 
expected to be accurate within ~ 15% at 7 mm. 

3. Observational results 

The continuum maps obtained at 0.87 mm, 1.3 mm, 
2.8 mm and 7 mm are shown in Fig.[T] Channel maps, inte- 
grated intensity and mean velocity maps of the ^^CO J=2- 
1, ^^CO J=l-0 and ^^CO J=3-2 transitions are shown in 
Fig.[2l[3]and[4]respectively. At the wavelength correspond- 
ing to the C^^O J=l-0 transition we did not detected any 
emission; the continuum integrated fluxes are reported in 
Tab.m 

^ http://cfa-www.harvard.edu/ cqi/mircook.html 
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A 


beam 




emission 


flux 


(mm) 


FWHM 


PA 


FWHM 


(mJy) 


0.87 


3". 14 X 2". 21 


20° 


3". 60 X 2". 64 


1910±20 


1.3 


1".95 X 0".42 


7° 


2". 01 X Q".83 


705±12 


2.8 


3". 3 X 0".94 


8° 


3". 5 X 1".4 


77.Q±2.2 


7 


1".71 X 0".81 


172° 


1".71 X 0".94 


4.5±0.5 



Table 1. Column 2 and 3 show the size and the position 
angle of the synthesized beam of the HD 163296 observa- 
tions performed with SMA (at 0.87 mm), PBI (at 1.3 mm 
and 2.7 mm) and VLA (at 7 mm). Column 4 show the 
FWHM size of the detected emission along the major and 
minor axis of the beam. Column 5 shows the integrated 
flux with the uncertainties due to the statistical errors. 



With simple physical assumptions, these observations 
allow us to determine some fundamental parameters of the 
star+disk system as the stellar and the disk masses, the 
contribution to the observed fluxes of free-free gas emis- 
sion, the wavelength dependence of the observed flux at 
millimeter wavelengths and the related dust grain opacity. 
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be obtained by comparing the observed velocities with the 
Kcplcrian law: 

v, = C-e-^l\ (1) 

where e is the angular distance from the central star and 
is the component of the disk rotational velocity along 
the line of sight. If the stellar mass M^, is in solar units, 
the stellar distance d in parsec and G is the disk incli- 
nation (0 = means pole-on disk), the constant C ~ 
30-\/A/*/(isin6' is the component along the line of sight 
of the disk rotational velocity (in km/sec) at e = 1". As 
shown in Fig.[5]the envelopes of both the ^^CO and ^^CO 
emissions are in agreement with C ~ 2.7 ± 0.4 km/sec. 
This value corresponds to a stellar mass of 2.0 ± 0.5 Mq 
for an inclination of 45° and d= 122 pc. Note finally that 
the disk rotation is clearly observed at least up to a dis- 
tance of about 4", corresponding to a minimum disk outer 
radius of about 500 AU. We will discuss in more detail the 
determination of the stellar mass and the disk outer ra- 
dius in Sec. [5l using a detailed model for the continuum 
and CO molecular emission. 



3.1. Disk morphology and apparent size 

The 12" X 12" continuum maps of HD 163296 are shown in 
Fig. [TJ At all wavelengths, the peak is coincident with the 
position of the optical star as measured from Hipparcos 
and, given the respective synthesized beams FWHM (see 
Tab. [1]), the emission is resolved and elongated approx- 
imately in the east-west direction. Approximating the 
source with a circularly symmetric geometrically thin disk 
and taking into account the beam shape, the observed as- 
pect ratio of the level contours implies an inclination of 
the disk plane from the line of sight of 45° ± 20° and 
a position angle of 120° ± 30°, in rough agreement with 
the values obtained by Mannings and Sargent (jl997p us- 
ing marginally resolved OVRO observations (58° and 126° 
respectively) . 

The emission is resolved at all the wavelengths in the 
East- West direction. The fainter contour levels are not 
centrally symmetric, showing excess emission in the east 
half of the image, better visible in the 0.87 and 1.3 mm 
maps. Both the disk size and the morphology will be dis- 
cussed in more detail in Sec. [H 

3.2. Disk kinematics 

As shown in Fig. [21 14] and [3l in all the detected molec- 
ular line transitions, the emission is resolved, showing a 
velocity pattern typical of an inclined rotating disk char- 
acterized by a position angle of about 130° (a more precise 
estimate of the position angle will be presented in Sec.[5|). 
The velocity-position diagram calculated along this direc- 
tion (see Fig. [5|) shows a well defined "butterfly shape" 
typical of Keplerian rotation. A first estimate of the mass 
of the central object and of the dimension of the disk can 



3.3. Free-free contribution and spectral index 

The continuum spatially integrated fluxes of HD 163296 
and the corresponding spectral energy distribution are 
given in Tab. [T] and shown in Fig. [6l respectively. In ad- 
dition to our measurements (full squares), observations at 
0.7, 1.3, 3.6 and 6.1 mm are from Natta et al. (2004 and 
references therein) while observations at 0.75, 0.8, 0.85, 
1.1 and 1.3 mm are from Mannings 111994^ . 

Assuming that the observed flux at 3.6 and 6.1 cm 
is dominated by free-free emission from a wind (with a 
spectral index of 0.6), the free-free contribution at 7 mm 
corresponds to 1.2 mJy (27% of the observed flux), while 
it is negligible at shorter wavelengths. 

After the subtraction of the free- free component, the 
millimeter spectral index Q!,„„i (-Fa oc ) calculated 

between 1 mm and 7 mm is 3.0 ± 0.1, slightly higher than 
the value amm — 2.6 ±0.2 obtained by Natta et al. (|2004p 
using VLA and OVRO fluxes only. 

3.4. Disk mass 

Assuming that the dust emission is optically thin at mil- 
limeter wavelengths and that the dust is isothermal, the 
measured flux can be used to estimate the product of the 
disk total mass M times the dust opacity ki, , through the 
relation: 



where d is the source distance and T is the temperature 
of the emitting dust. At millimeter wavelengths, the dust 
opacity can by parametrized by a power law 

fc = 0.01 • (A/1.3mm)-'3cm2 g-i, (3) 
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Fig. 1. Continuum maps of HD 163296 at 0.87, 1.3, 2.8 and 7 mm, starting form the left. In order to highlight the 
extended morphology of the disk, the first contour level In the 0.87 mm map corresponds to 30 mJy {3a), the second 
to 6tT while the inner contour levels are spaced by lOtr. At longer wavelength the contour levels are all spaced by 3<t, 
corresponding to 12 mJy at 1.3 mm, 3.3 mJy at 2.8 mm and 0.75 mJy at 7 mm. The small boxes show the relative 
synthesized beams. The integrated fluxes, the beam dimensions and orientations are summarized in Tab. [1] 



where the normalization at A = 1.3 mm assumes a 
dust/gas mass ratio of 0.01 (Beckwith et al.. fTMOl fTMTl) . 
Taking a characteristic temperature for the dust in the 
outer disk of an early A star of 30 K (Natta et al. I2000p ,we 
derive a total gas+dust mass of 0.12 Mq for the circum- 
stellar material in HD 163296. This rough mass determi- 
nation, that does not take into account the presence of 
optically thick emission from the inner part of the disk, 
nor of temperature gradients, will be discussed in more 
detail in Sec. 15.21 

With these simplifying assumptions, the slope of the 
dust opacity is/3 = a — 2^1, where a is the spectral 
index obtained in the previous section. 

4. Disk Models 

The disk parameters derived in Sec. [3] under a num- 
ber of very crude assumptions can only provide order- 
of-magnitude estimates. A more quantitative analysis re- 
quires to compare the observations to more sophisticated 
model predictions. We chose to perform the comparison 
using the observed visibilities (rather than reconstructed 
images) , following a similar approach as the one discussed 
by Dutrey et al. ((2006;) . 

This method requires to decide "a priori" which fam- 
ily of models is likely to describe the observed object. 
In view of the results described in Sec. [3l we model the 
millimeter emission of HD 163296 (continuum and CO 
lines) as coming from a circumstellar disk. We assume 
that the disk is heated by the stellar radiation only, and 
that any viscous contribution can be neglected. This is 
very likely a good approximation, given the relatively low 
accretion rate measured in HD 163296 (~ 10~^ Mo/yr; 
Garcia Lopez et al.. 12006 1). 

The disk structure and emission is computed using the 
2-layer approximation of Chiang and Goldreich (jl997p . as 
developed by Dullemond et al. (|200ip . Similar models have 



been used in Testi et al. (2003) and Natta et al. (|2004l) . 
to analyze the (sub)millimeter emission of a number of 
Herbig Ae stars. We refer to these papers for a more de- 
tailed description. 

The inner disk is characterized by the presence of a 
puffed-up inner rim located at the dust evaporation radius 
(Isella and Natta. l2005p . The combined rim-|-disk flux is 
calculated taking into account the shadow that the rim 
casts over the outer disk. 

The disk geometry can be fully flared, as in hydrostatic 
equilibrium disks where gas and dust are well mixed. If 
dust growth and settling are important, the flaring an- 
gle decreases: the resulting radial temperature profile is 
flatter, affecting the SED in the mid and far infrared 
(D'Alessio et al. 120061 and references therein). As shown 
in Fig. [51 the SED of HD 163296 is weh reproduced by a 
fully flared disk model. The fit of the SED becomes very 
poor if the flaring angle is reduced by more than a factor 
of two, which corresponds to a small variation on temper- 
ature radial profile. 

Once the stellar properties are known, the disk struc- 
ture is completely characterized by the following parame- 
ters: the disk mass (M^), the disk outer radius (Rout), the 
dependence of the surface density on radius (S cx r~P) 
and the properties of dust on the disk surface and mid- 
plane. In addition, the observed emission depends on the 
orientation of the disk with respect to the observer, which 
is characterized by the inclination t of the disk with re- 
spect to the line of sight (t = for face-on disks) and the 
position angle PA. 

4.1. Continuum emission 

The continuum emission at millimeter and sub- millimeter 
wavelengths is computed by ray integration as in 
Dullemond et al. (|200ip . We describe the midplane dust 
opacity at long wavelengths as a power law of index /3, 
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Fig. 2. Velocity channel maps of the ^^CO 3=2-1 emission. The velocity resolution is 0.25 km/sec. The LSR velocity 
is indicated in the upper left corner of each panel. The angular resolution (synthesized beam) , indicated in the small 
boxes, is 1.95" x 0.42" at PA 7°; the contour spacing is 0.23 Jy/beam corresponding to 3a. The last two panels show 
the integrated intensity (contour levels spaced by 0.4 Jy/beam km s~^) and the velocity field (contour levels from 
3 km s~^ to 9 km s~^ spaced by 0.5 km s~^) respectively. 



with (3 a free parameter, as in Eq.[3] At short wavelengths 
and in the disk surface we adopt the dust opacity of as- 
tronomical silicates (Weingartner and Draine. r200ip . The 
inner disk radius is the dust sublimation radius, as in the 
rim models of Isella et al. (|2006p for large (~ lyum) grains, 
and it is equal to 0.45 AU; the dust on the disk surface is 
as in Natta et al. (|2004p . Neither of these two quantities 
is relevant for the following analysis. 



4.2. CO emission 

The observed CO emission originates in the outer layers 
of the disk, at heights that depend on the optical depth of 
the specific transition. Once the disk structure is specified, 
as described above, one needs to compute, at each radius, 
the gas temperature profile in the vertical direction. This 
is a complex problem, whose results depend on a num- 
ber of not well known properties, among them the X-ray 
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Fig. 3. Velocity channel maps of the ^^CO J=l-0 emission. The velocity resolution is 0.25 km/sec. The LSR velocity 
(km s^^) is indicated in the upper left corner of each panel. The angular resolution (synthesized beam), indicated in 
the small boxes, is 3.3" x 0.94" at PA 8°; the contour spacing is 0.09 Jy/beam corresponding to 3a. The last two 
panels show the integrated intensity (contour levels spaced by 0.12 Jy/beam) and the velocity field (contour levels 
from 3 km s^^ to 9 km s~^ spaced by 0.5 km s~^) respectively. 



field and the role of very small grains in heating the gas 
(e.g., DuUemond et al. I2006p . Therefore, we use a para- 
metric description assuming that for each CO transition 
the excitation temperature is the same at all z and can be 
described as a power-law of r in the form: 



(4) 



The assumption of a constant excitation temperature 
along the vertical direction is correct for optically thick 



lines, and/or if the velocity gradient along the line of sight 
is very large. In both cases, the contribution at each wave- 
length from the line of sight under consideration comes 
from a small region only, where the optical depth in the 
line is of order unity. For any CO line, the height to which 
Tiine{f) refers is therefore different, depending on the den- 
sity and velocity structure, as well as on the disk inclina- 
tion angle. The values oiTuneiro) and q for each CO tran- 
sition are free parameters. A similar procedure has been 
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Fig. 4. Velocity channel maps of the "'^^CO J=3-2 emission. The velocity resolution is 0.18 km/sec. The LSR velocity 
(km s~^) is indicated in the upper left corner of each panel. The angular resolution (synthesized beam), indicated in 
the small boxes, is 3.4" x 2.21" at PA 20°; the contour spacing is 1.5 Jy/beam corresponding to Scr. The last two panels 
show the integrated intensity (contour levels spaced by 5 Jy/beam km s^^) and the velocity field (contour levels from 
3 km s^^ to 9 km s~^ spaced by 0.5 km s""'^) respectively. 



used by Dutrey et al. ([TOM|) and Dartois et al. (PUUH]) in 
their study of the CO emission of T Tauri disks. 



In the analysis of CO observations we have as- 
sumed 12CO/H2 = 7.0 • 10-^ 13CO/H2 = 1.0-10-6 and 
C^^0/H2 =1.3T0~^, which correspond to an isotopic ra- 
tio i2co£3co = 70, i2co/C^ = 550 (Beckwith and 
Sargent. 119931 Dutrey et al.. 119961 and reference therein) 



As described in detail in the Appendix, we have de- 
veloped a code which computes, for each CO transition, 
the line intensity and profile as function of the disk pa- 
rameters, namely the inclination and PA, the density and 
temperature profile, the disk outer radius. In addition, we 
assume that the disk is in Keplerian rotation around the 
central star and vary the stellar mass independently for 
each line. 



8 Isella A. et al.: Millimeter imaging of HD 163296 



Parameters 


Continuum 


-CO J=2-l 


-co J=3-2 


'^CO J=l-0 


(M©) 


2.6" 


^•"-0.5 


2.4 ± 0.8 


2.6 ± 0.6 


PA 


120° ± 20° 


128°± 5° 


130° ± 13° 


130°± 8° 


Incl 


40°± 12° 


45° ± 5° 


45° ± 10° 


50° ± 8° 


r> /ATTN 

Rout (AU) 


200 ± 15 


550 ± 50 


550 ± 100 


500 ± 80 


SioAu (g/cm^) 


46 ± 4 


90 ± 70 


90" 


41- 


P 


0.8 ± 0.1 


0.6t°:? 


0.6" 


1.0 ± 0.5 


Tco,woAu (K) 






60 ± 20 


30± 10 








0.8± 0.4 


0.8± 0.4 


P 


1.0 ± 0.1 


1.0" 


1.0" 


1.0" 




2.6 


1.17 


1.05 


1.15 



Table 2. Parameters of the disk structure relative to the best fit models for the continuum and the CO emissions as 
described in Sec. [3 For each parameters uncertainties are given at la level. ° Fixed parameter. 



In the following, we assume a gas turbulence velocity 
vturb = km/sec. Note that for Vturb in the range observed 
in other TTS and HAe (0.07-0.38 km/sec from Pietu et 
al. [20071 [20051 [20081 Dutrey et al. [T9M1 [T9961 Dartois et 
al. 120031 Simon et al. I200ip the results does not change. 
This is due to the fact that, for the inclination of the 
HD 163296 disk, the CO lines width is dominated by the 
differential disk rotation. 

5. Results 

5.1. Method of analysis 

The observations of HD 163296 have been analyzed by 
comparing the observed and the model predicted complex 
visibilities. For each set of data, the best fit model has been 
obtained minimizing the = ^i{{Re.mod,i - Reobs,if + 
{Inimod,i — I'mobs,iY) ' Wi, whcrc Rex and Irax are the 
real and the imaginary part of the complex visibility mea- 
sured {obs) and predicted {mod) in the point i of the uv- 
plane, and Wi is the weight of each measure. For each CO 
transition, the has been computed by the simultaneous 
fitting of 14 velocity channels chosen in order to optimize 
the sampling of the line profile. The minimization has 
been performed exploring a wide region of the space of 
model parameters. For each parameter, the Icr uncertain- 
ties are estimated as Xio- = Xm + V^^, where n is the 
number of degrees of freedom and Xm i^ the value of 
best fit model. 

5.2. Continuum emission 

The continuum data at the four observed wavelengths 
(0.87, 1.3, 2.7, 7 mm) have been analyzed independently 
to estimate the best disk parameters and their uncertain- 
ties. We have varied the two parameters which define the 
position of the disk on the plane of the sky, namely the 
inclination and position angle, as well as the three physi- 
cal parameters which affect the continuum dust emission, 
namely the disk outer radius Rout , the slope p of the sur- 



face density profile and the product Tiiqau ^ k\. Since 
the continuum emission has a very weak dependence on 
the mass of the central star, we fixed M^,=2.6 Mq (see 
Sec.lOl). 

The values of the parameters obtained at different 
wavelengths are the same within the uncertainties. In col- 
umn 1 of Tab [2] we show the values of the resulting best fit 
model obtained by combining the independent results at 
the four different wavelengths. The parameters constraint 
is dominated by the 1.3 mm data, which have the best 
sensitivity and resolution. 

The surface density radial profile has a slope p = 
0.81±0.01 and the disk outer radius is Rout = 200±15 AU. 
The very good constraints on both parameters is due to 
the favorable orientation of the beam at 1.3 mm with the 
maximum resolution (0.42") in almost the same direction 
of the major axis of the disk. Note that the formal error 
on p is extremely low. On the other hand, small varia- 
tions on the disk flaring (that allow to fit the SED) lead 
to correspondingly small variation of the disk radial tem- 
perature profile and on the value of p. Taking into account 
this fact, ±0.1 is a more reasonable uncertainty for p. The 
formal uncertainties on the other parameters are much 
larger than those introduced by small variations of the 
flaring angle. 

The four values of T^ioau x k\ are used to constrain the 
slope (3 of the dust opacity law (see Eq. [3|) between 0.87 
and 7 mm, which turns out to be 1.0 ± 0.1. This confirms 
the presence of large grains in the HD 163296 circunistellar 
disk (Natta et al. l2004p . Note that the optically thick disk 
region at 1.3 mm has size of '--^ 6 AU and contributes 
only 4% of the observed flux, so that 13 is very similar to 
that derived in Sec. 13.31 The observations constrain only 
the dependence of the dust opacity on wavelength, not 
its absolute value. The value of Sioa;7 given in Tab. O is 
obtained from Sioac/ x fci.snim assuming /ci.smni = 0.01 
cm^/g. 

The residuals for each wavelength are shown in Fig. [71 
They have been reconstructed from the residuals in the uv- 
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R.A. 

Fig. 7. Maps of the residuals relative to the best fit model for the continuum emission (see Fig. [1] for the observations). 
The contour level are spaced by 3a corresponding to 30 mJy at 0.87 mm, 12 mJy at 1.3 mm, 3.3 mJy at 2.8 mm and 
0.75 mJy at 7 mm. The small boxes show the relative synthesized beams. 



plane with the same procedure used to obtain the images 
in Fig. [TJ Residual contours are generally lower than 3(t 
with the exception of the 0.87 and 1.3 mm maps, where 
a flux asymmetry in the east half of the map (see also 
Sec. 13. ip is clearly visible. This structure, not detected at 
longer wavelengths, requires more resolved observations to 
be investigated in detail. 

5.3. CO emission 

The analysis of the CO emission has been carried out sepa- 
rately for the different CO transitions and the correspond- 
ing best fit parameters are given in Tab. [2l 

The observed "'^^CO J=2-l emission is well fitted by 
a Keplerian disk orbiting a central star with mass Mi, = 
2.6lg;| Mq. Fig. [5] shows the position-velocity residuals 
obtained subtracting the best fit model from the observed 
uv-table: no evidence of non-Keplerian rotation or stellar 
outflow is detected, within the actual instrumental sensi- 
tivity. Both the position angle and the inclination are in 
agreement with the values obtained from the continuum. 
Since the line is optically thick, the constraint on the gas 
surface density is poor. The inferred outer radius of the 
disk is 550 ± 50 AU, more than two times larger than the 
value obtained from the continuum and similar to the re- 
sult obtained by Thi et al. (|2004p from the model fitting 
of single dish line profiles. The radial temperature profile 
of the i^CO has a slope q = O.^toi and a value of 40^5 K 
at 100 AU. 

Since the ^^CO J=3-2 line is optically thick and given 
the relatively low resolution of the SMA observations, the 
measurements do not constrain the CO radial density pro- 
file. In this case we fix both p and "Siqau equal to the 
values obtained for the ^^CO J=2~l and vary the other 
parameters. The results for the inclination, position angle, 
disk outer radius and CO temperature are in good agree- 
ment with the values obtained for the ^^CO J=2-l, with 
larger uncertainties due to the lower spatial resolution. 

Finally, the model fit to the ^^CO J=l-0 line gives re- 
sults consistent with those obtained from the ^^CO lines 



with the exception of the value of Eioa[/ which is signifi- 
cantly smaller. This discrepancy may be due to a depletion 
of the ^"^CO and will be discussed in Sec. 16.51 The radial 
temperature profile of the ^'^CO has a slope q — 0.8 ± 0.4 
and a value of 30 ± 10 K at 100 AU. We will comment 
on the gas physical conditions in Sec. 16.51 In general, the 
parameter constraints obtained from the ^^CO are not so 
good as those obtained for the ^^CO J=2-l, as expected 
given the lower resolution of the observations. 

6. Discussion 

6.1. Disk orientation 

The disk orientation is well determined by our set of mea- 
surements: it has it has a moderate inclination with re- 
spect to the line of sight (the mean of the inclination val- 
ues is 46° ± 4°) with a position angle of 128° ± 4°. 

6.2. Disk kinematics and stellar mass 

All the observations are consistent with the emission of 
a circumstellar disk in Keplerian rotation around a star 
of 2.6 M0, assuming the Hipparcos distance of 122j;]^3 pc 
(van den Ancker et al. ll998p . Within the error, the stellar 
mass is in agreement with the value of 2.3 Mq (Natta et 
al., I2004P obtained from the location of the star on the 
HR diagram, using Palla and Stabler (|1993p evolutionary 
tracks; the corresponding stellar age is of about 5 Myr. 

6.3. Disk outer radius 

The model fitting (Tab. [2]) shows that the value of the 
disk outer radius inferred from the continuum dust emis- 
sion (200 ± 15 AU) is almost three times smaller than the 
value obtained from the CO analysis (where 540 ± 40 AU 
is the mean value). Since our method takes into account 
the sensitivity limits of the different interferometric ob- 
servations, this discrepancy can not be explained by the 
fact that the disk outer regions have a continuum surface 
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Fig. 5. Velocity-position plots along the plane of the disk 
for the 12 CO J=2-l (upper panel), the ^^CO 3=3-2 (mid- 
dle panel) and the ^^CO J=l-0 (lower panel) transitions. 
The angular offset is measured with respect to the phase 
center of the observations corresponding to the position of 
the central star. The contour levels are spaced by 2cr corre- 
sponding to 0.14 Jy/beam, 1 Jy/beam and 0.06 Jy/beam 
respectively. The cross in the lower left of each panel gives 
the angular and spectral resolution of the correspond- 
ing map. The thick solid lines marks the border where 
emission is expected for a Keplerian disk inclined by 45° 
and rotating around a 2.0 M0 point source; the exter- 
nal and internal dashed lines correspond to stellar masses 
of 2.5 M0 and 1.5 Mq respectively . The horizontal and 
vertical straight dashed lines mark the systemic velocity 
(5.8 km/sec) and the position of the continuum peak. 



brightness below the sensitivity limit. If we extend the 
disk model that fit the continuum to the outer radius of 
the CO, we predict a continuum emission between 200 AU 
and 550 AU which would be easily detected in our observa- 
tions (see Fig. |9]). A similar difference in the dust and gas 
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Fig. 6. HD 163296 spectral energy distribution. 
Data from literature are shown with empty squares 
(Mannings [UMl Natta et al. [200l de Winter et aL [200T|) 
and dots (ISO-SW). Our new measurements (see TabHJ 
are shown with full squares. The observed fluxes at 3.6 cm 
and 6.1 cm have been used to calculate the free- free con- 
tribution (dashed line) to the observed flux assuming a 
free-free spectral index of 0.6. Free-free subtracted fluxes 
are shown at 7 mm while at shorter wavelength the free- 
free contribution is negligible. The solid line shows the 
prediction of our disk model as discussed in Sec. 15.21 The 
resulting spectral index amm between 0.87 and 7 mm 
(Fa oc A^"""") is 3.0 ±0.1. 



outer radii has been also reported by Pietu et al. (|2005p 
for the Herbig Ae star AB Aur. 

To reconcile CO and dust observations, it is necessary 
to introduce a sharp drop in the continuum emission of 
a factor > 30 at a radius of about 200 AU. With such 
a drop, the millimeter fluxes at larger r will be below 
the sensitivity of our observations and will be lost in the 
observational noise. What can be the origin of such a drop? 

In the optically thin regime, which is appropriate for 
the outer regions of the disk, the continuum flux emitted 
at distance r from the star depends on the mass surface 
density S(r), the dust/gas ratio 11 (r), the dust opacity 
ki,{r) and the midplane dust temperature T(r) through 
the relation 

F,(r)cxS(r)-n(r)-fc,(r)-r(r). (5) 

It is possible that a different disk geometry (e.g., a 
lower flaring angle) gives a midplane dust temperature 
lower than the values predicted by our disk model. A lower 
temperature limit, of ~10 K, is however imposed by the 
equilibrium with the interstellar radiation fleld. Given that 
our fully flared disk model predicts temperatures of 20- 
30 K in the millimeter emitting regions, one can reduce 
Fi, by a factor 3 at most. 

A second possibility to explain the observed flux de- 
pletion is that the dust opacity k,y{r) at distance larger 
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Offset (arcsec) 

Fig. 8. Comparison between the observed and the model 
predicted ^^CO J=2-l emission. The upper panel shows 
the position-velocity diagram for the ^^CO J=2-l transi- 
tion (as in Fig. [5]). The lower panel shows the residuals rel- 
ative to the best fit model parameters reported in Tab. [2] 
The contour levels are at 2a as in Fig. [51 



than 200 AU is much lower because most of the grains 
have grown into very large bodies. 

Finally, a third possibility is that the dust density, i.e. 
I](r) • n(r), is in some way depleted at large distances. 
The values of the gas density E(r) obtained by the CO 
lines analysis (see Tab. induce however to exclude the 
existence of the strong discontinuity in the gas radial den- 
sity profile required to explain a flux depletion factor > 30 
since a drop of a factor ~10 in the ^^CO distribution at 
T > 200 AU will already produce small but observable ef- 
fects on the line emission. On the other hand, it is possible 
that the ratio dust/gas n(r) decreases very rapidly in the 
outer disk either due to the formation and fast migration 
of meter size bodies or to the presence of planetesimals 
or planets which create a gap in the dust distribution. It 
is worth to note that the dynamical perturbation induced 
by such large bodies on the surrounding material (see the 
review of Papaloizou et al. I2007P may also account for the 
asymmetric dust density distribution detected in the con- 
tinuum maps (Fig. [7]) . 

The presence of a giant planet, or a brown dwarf, or- 
biting in the outer disk of HD 163296 has also be sug- 
gested by Grady et al. (|2000p . in order to explain the dark 
line observed in the scattered light HST images between 
300 AU and 350 AU from the central star. While planets 
are invoked to explain the large inner gaps observed in the 
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Fig. 9. Model predicted 1.3 mm continuum emission ob- 
tained extrapolating the dust surface density inferred in- 
side the radius of 200 AU to the CO outer radius of 
550 AU. The contours are spaced by 3(7 corresponding 
to 0.12 mJy. See Fig. [T] for the comparison with the ob- 
servations. 

dust distribution in the so called "transitional disks" (i.e., 
Calvet et al, [M}5|) . HD 163296 will be, if confirmed by 
future observations, the first case in which a sub-stellar 
mass companion is found in the outer disk of a pre-main 
sequence star. 

6.4. Disk surface density and mass 

Both the dust continuum emission and the CO lines in- 
dicate a rather shallow surface density profile, E oc r~°-* 
corresponding to a gas -I- dust mass of 0.05 ± 0.01 M0 if 
the disk outer radius is 200 AU. However, if, as we believe, 
the gaseous disk extends with the same surface density 
profile all way to 540 AU the inferred disk mass is much 
larger, 0.17 Mq. 

An additional uncertainty on the disk mass comes from 
the uncertainty on the dust opacity, as discussed in Natta 
et al. (|2007|) . Although a disk mass of 0.17 Mq(6% M^) 
is probably a lower limit, the observed Keplerian rota- 
tion pattern indicates that the disk should not be much 
more massive. In this case one expects deviations from 
Keplerian rotation (Lodato and Berlin, I2003P as observed 
in massive disks (Cesaroni et al.. I2005|) . 

The example of HD 163296 illustrates how the simple 
disk mass derivation from the observed millimeter fiux (see 
Sec 13.41) can be wrong for a variety of reasons, even when 
the emission is optically thin. In the case of HD 163296, 
only 30% of the total mass (i.e, the fraction inside 200 
AU) contribute to the continuum millimeter fiux. On the 
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Fig. 10. Comparison between the dust and the CO tem- 
perature. As labelled in the figure, the thin-solid line cor- 
respond to the dust temperature in the optically thin disk 
surface, while the thick solid line corresponds to the dust 
temperature in the disk interior. The CO temperature is 
indicated with the shaded region, taking into considera- 
tion the uncertainties. The model parameters are those 
reported in column 3 of Tab. [5] 



other hand, the disk is hotter than the 30 K assumed in 
Eq. [21 The two effects cancel in part so that the "simple" 
and "correct" values differ by about 50% (for the same 
fci.amm normalization). 

6.5. CO vs dust temperature 

At the actual gas density both the ^^CO and ^■^CO 
lines are optically thick and lead to similar CO exci- 
tation temperatures expressed by the relation T{r) ~ 
40K-(r/100AU)~°-^. Fig. [10] shows the comparison be- 
tween this CO temperature (shaded region) and the dust 
temperature corresponding to the same disk model. The 
disk interior temperature Ti (thick-solid line) and the disk 
surface temperature Tg (thin-solid line), calculated using 
our two-layer disk model, correspond respectively to disk 
regions where r^, >> 3 and r^, ~ 1 (DuUemond et al. l2002p : 
Ti, is the optical depth for the stellar radiation calculated 
along the radial direction. The equivalence between the 
CO and the dust temperature can thus be interpreted as 
the evidence that the CO, and in particular the ^^CO, 
lines are emitted by the gas present at such values of the 
optical depth, e.g. well under the disk surface. Is this re- 
sult compatible with the high optical depth of the ^^CO 
lines? 

Assuming that H2 and CO are well mixed all over the 
vertical extension of the disk and assuming also the stan- 
dard ratio l2cO/H2=7.5•10-^ the ^^CO lines should be 
emitted by disk regions well above the disk surface, de- 
fined as the height where = 1. In this situation, we 
expect a CO temperature close to the disk surface temper- 
ature, or even higher if the gas and the dust are decoupled 



(Jonkheid et al., 2006). On the other hand CO molecules 
on the superficial layers of the disk are expected to be 
photo-dissociated by the UV photons emitted by the hot 
central star. In particular, for an Herbig Ae star similar 
to HD 163296, Jonhkeid et al. (2006) show that the CO 
photo-dissociation occurs in the disk regions where the 
optical depth to the stellar UV radiation field is tuv ^ 3 
(see Fig 5.6 therein). Therefore, the CO and disk inte- 
rior temperature are expected to be similar if rj/y < r*, 
i.e., if kjjv ^ , where is the dust opacity averaged 
on the stellar spectrum. We think that this is the case in 
HD 163296. 

The dust opacity is a complex quantity that depends 
on the dust chemical composition, structure and on the 
grain size distribution. The condition kjjv ^ is gen- 
erally satisfied by large grains. In the case of HD 163296 
one needs, for example, compact silicate grains larger than 
0.1/im or porous grains of silicate and carbonaceous mate- 
rials larger than few microns (see Fig. fTTj) . Cooler stars re- 
quire larger grains, unless significant UV excess is present. 

If CO survives only in the disk interior, we expect that 
lines with different optical depth will have similar temper- 
atures. On the contrary, if CO is present in the disk surface 
layers, lines with higher optical depth should have higher 
excitation temperatures. 

The few existing observations indicate that all cases 
occur. Pietu et al. (2007) have recently reported that for 
the HAe stars MWC 480 and AB Aur and the TTS DM 
Tau, the ^^CO and ^'^CO lines are characterized by differ- 
ent temperatures, while in the case of the TTS LkCa 15 
the CO temperatures are all similar. We suggest that these 
differences are due to different dust properties (i.e. differ- 
ent composition, grain growth and settling) and stellar 
spectra. 

6.6. CO isotopic abundances 

The values of SioAt/ reported in Tab. [2] indicate that the 
^^CO J=l-0 emission requires a ratio ^'^C0/H2 ^ 10^^, 
about a factor 10 lower than what found in interstellar 
clouds. From the non-detection of the C^'^O J=l-0 line we 
estimate that the C^^O should be depleted of a factor > 60 
with respect to the typical ratio ^^CO/C^^O = 550. We 
note that a similar trend of depletion from ^^CO to "'^^CO 
and C^^O has been found by Dutrey et al. (fT994l [T9961) in 
a number of TTS. On the other hand an over-abundance of 
"CO has been found in DM Tau, MWC 480 and LkCa 15 
(Pietu et al. 2007) while for DM Tau, Dartois et al. (2003) 
find that all the CO isotopomers (^^CO, ^^CO and C^^O) 
are depleted by the same factor ~10. The observational 
picture is still very uncertain and deserves further inves- 
tigations. 

In the case of HD 163296, the gas temperature always 
higher that 20 K roles out the condensation of CO onto 
dust grains as possible cause of the CO isotopomers de- 
pletion. Alternatively, Jonkheid et al. (|2006p have recently 
pointed out that an higher ratio between the ^^CO and its 
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Fig. 11. Ratio between the Planck mean opacity at the 
stellar temperature and the UV opacity kuv for differ- 
ent grain sizes and composition. The upper panel shows 
the opacity ratio for compact silicate grains; the lower 
panel shows the ratio for porous grains composed by 
Olivine (17%), carbonates (53%) and vacuum (30%). The 
solid lines are relative to an HAe star of ~ lOOOOK, 
the long-dashed lines refer to a TTS with = 4000 K 
and the short-dashed lines to the same TTS with the 
addition of the UV excess radiation calculated so that 
Lijv/Li, = 0.05. This latter case is representative for the 
emission released by the gas accreting on the stellar sur- 
face. 



isotopomers may be the consequence of grain growth and 
settling. This latter hypothesis is more convincing given 
the evidence of grain growth in the HD 163296 disk dis- 
cussed in this work (see Sec. I5.2landl6.5ll. 



7. Summary 

This paper presents new observations of the disk of 
HD 163296 in the dust continuum from 0.87 to 7 mm, 
i^CO (J==2-l and J=3-2) and ^^CO (J=l-0) lines. The 
disk is resolved in all lines and continuum. 

We have compared the observations to the predictions 
of self-consistent disk models. We find that the disk, as 
seen in CO lines, is very large [R = 540 AU), with a 
Keplerian rotation pattern consistent with a central mass 
of 2.6 M0. Within the observational errors, there is no 
evidence of non-Keplerian motions and/or significant tur- 
bulent broadening. We obtain a disk inclination of 46°, 
significantly lower than the value of ~60° found in litera- 
ture. 

The dust opacity has a power law dependence on wave- 
length K oc with /3 = 1.0 ± 0.1 in the interval 0.87-7 
mm. This value is similar to what has been measured in a 
number of spatially resolved disks of HAe and TTS (e.g., 
Natta et al. 120041 Rodmann et al. l2006p . and is very likely 
an indication that the bulk of the solid material in these 
disks has coagulated into very large bodies, of millimeter 
and centimeter size (Natta et al. I2007P . Within the accu- 



racy of our data we do not confirm the possible variation 
of (3 with r discussed in Natta et al. ()2007p . 

The continuum observations constrain the surface den- 
sity profile (S cx r°-^) for r < 200 AU. At larger radii, the 
continuum emission drops with respect to the model pre- 
dictions by a factor 30 at least (at 1.3 mm). We argue 
that this may be due to the clearing of a very large gap 
by dynamical perturbations from a low mass companion 
or to the inward migration of the large bodies that may 
form in the outer disk. 

The temperatures derived for CO lines of different op- 
tical depth are similar and equal to the dust temperature 
in the disk interior (r^, >> 1). This requires that the dust 
opacity in the UV (which controls the CO dissociation) 
and in the wavelength range where the stellar radiation 
peaks are similar, as expected if grains have grown to mi- 
cron size. Differences in the temperature of the CO lines 
can only be expected if CO is present in the disk surface 
layers. A variety of situation can occur depending on the 
exact dust composition and sizes and on the stellar radi- 
ation field. 

A comparison of the disk properties derived from the 
dust continuum and the CO lines shows that the ^^CO 
J=l-0 emission is consistent with a ratio of ^'^CO/H and 
C^^O/H about a factor 10 and 60 lower than what is found 
in interstellar clouds. 



8. Conclusions 

The results discussed in this work can be interpreted as 
clues of the evolution occurring in the HD 163296 system. 
The presence of large grains in the disk midplane, the 
equivalence between the dust and CO temperature, the 
drop in the continuum dust emission further out 200 AU, 
its detected asymmetry and, maybe, the CO isotopomers 
depletion, all support the idea that the circumstellar disk 
is probably harboring the formation of large bodies, being 
in between a Class II pre-main sequence disk and an older 
debris disk. 

On the other hand, it is important to underline that 
HD 163296 strongly differs from the so called transitional 
disks (Calvet et al., 2005) characterized by a clearing of the 
inner disk supposed to be originated by the presence of a 
giant planet. In the case of HD 163296 is indeed the outer 
disk that appear dust depleted. This conclusion opens 
a number of theoretical problem about how dust grains 
can growth in the low density outer disk environment. In 
this respect the recently improved PdB array, the new 
CARMA array and ALMA will bring in the next future 
to real observational breakthrough in the comprehension 
of circumstellar disk evolution and planetary formation. 

Acknowledgements. The authors acknowledge partial support 
for this project by MIUR PRIN grant 2003/027003-001. A.I. 
acknowledge Riccardo Cesaroni, Malcolm Walmsley, Giuseppe 
Bertin, Ewine van Dishoek and Anneila Sargent for their help 
and useful suggestions. 



14 



Isella A. et al.: Millimeter imaging of HD 163296 



Appendix A: CO emission model 

In order to analyze the CO emission we developed a nu- 
merical code that solve the general formulation of the ra- 
diation transfer equation along each direction between the 
observer and the emitting source. If s is the linear coordi- 
nate along the line of sight, increasing from the observer 
(s = 0) towards the source, the observed emission in each 
direction is given by the relation 



(A.l) 



where the optical depth is defined in each point s through 



Ku{s')ds' 



(A.2) 



and is the absorbing coefficient of the interstellar 
medium. Given the high gas densities in the pre-main se- 
quence circumstellar disks, we can assume that all the 
CO levels corresponding to the rotational transitions un- 
der investigation are thermalized. In this case the source 
function can be approximated by the Planck function 



2hi/3 



1 



exp(hj./kTco(s)) - 1 



, (A.3) 



depending only on the local temperature of the gas Tco 
(c, h and k are respectively the light speed, the Planck 
and Boltzman constant). 

The absorbing coefficient of the circumstellar medium 
is due both to gaseous CO and dust: Ky{s) = K^^{s) + 
K^{s). For the dust K^{s) = p{s) ■ k^,, where p{s) is the 
local density of the circumstellar material (gas-|-dust) and 
kv is the dust absorbing coefficient for total mass unit 
given in Eq. [3l The CO absorbing coefficient is given by 
the relation 



where ni{s) is the total number of CO molecules at the 
lower level I of the transition and (Ji/{s) is the CO absorb- 
ing cross section. 

Calling mo and Xco the mean molecular weight of the 
gas and the fraction of CO present in the gas respectively, 
the number of molecules n;(s) is given by the Boltzman 
equation 



ni{s) = xco — TT^' 

mo Z{Tco{s)) 



(A.5) 



where gi = 2Z -|- 1 is the statistical weight of the lower level 
I of the transition, Ei = (l/2)Z(/-|-l)kri is the level energy, 
Ti is the temperature equivalent to the transition energy 
and Z{Tcois)) is the partition function at the gas tem- 
perature Tco(s). Following Beckwith and Sargent p9931 
and references therein), the absorbing cross section a^{s) 
can be expressed in term of the integrated cross section of 
the transition (Tq through the relation 



a.(s)-f7o-0.(s)-(l-e-h'^/kTco(s))^ 




Fig. A.l. Schematic representation of the frame of ref- 
erence adopted to calculate the CO emission arising from 
a rotating disk. The disk midplane and the observer lie 
respectively on the {x, y) and {y, z) planes; 8 is the disk 
inclination; d is the distance between the observer and 
the central star; s is the linear coordinate along the line 
of sight increasing from the observer (s = 0) towards the 
emitting source. Assuming that the material within the 
disk is subject to the Keplerian rotation around the cen- 
tral star, we call the velocity of a mass element at 
distance r and Vk the projection of ti^ along the line of 
sight. 



where 4>^{s) is the intrinsic line profile 

c 1 / A^v 

I/O AVV^ 

with 



exp 



A2y 



(A.4) AV 



and 



l2kTco{s) 
mco 



S^^kTi {I + 1)2 
h^c 



"turb 



(A.7) 



(A.8) 



(A.9) 



21 1 

In the previous equations, mco is the CO molecular 
weight, Vf) is the rest frequency of the molecular tran- 
sition, Au is the difference between the velocity Vobs = 
{c/vq){v — i^o), corresponding to the frequency v, and the 
component of the gas velocity along the line of sight, i'fe(s), 
ji is the dipole moment of the CO molecule. Note that 
writing Eq. IA.7[ we assumed that the intrinsic line width 
depends both on the thermal velocity dispersion in the gas 
and on turbulent velocity vturb- 

As shown in Fig. [5] and described in Sec. 13.21 the ob- 
served velocity patterns in the CO transitions are in good 
agreement with the Keplerian rotation of the disk. We can 
thus assume that gas moves on circular orbits around the 
star characterized by a tangential velocity 



(A.6) 



CM. 



(A.IO) 
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where r is the radius of the orbit and is the stellar 
mass. To calculate the velocity i'fc(s), component of the 
tangential velocity v^{r) along the line of sight, it is use- 
ful to define a coordinates system centered on the star as 
shown in Fig. lA.ll the (x, y) plane corresponds to the disk 
midplane; the observer lies in the (y, z) plane and its po- 
sition it is defined by the inclination 6 and the distance d 
from the star; each point of the circumstellar spa ce can be 
defined through the cylindrical coordinates r — \J + J/^, 
= arctan {y /x) and z. Since r, we can write 



Vk{s) = Vcf,{r) ■ cosip ■ sin^ 



(A.ll) 



Finally, the velocity Vk{s) can be calculated for each di- 
rection, knowing the geometrical transformations between 
the coordinate s along the line of sight and the cylindrical 
coordinates r and 0. 

In order to solve the described set of equations, we thus 
need an expression for the circumstellar mass density p{s) 
and the temperature of the emitting gas Tco{s)- In both 
cases we can assume the cylindrical symmetry and write 
p{s) = p{r, z) and Tco{s) = Tco{r)- 

For the emitting gas temperature, we choose the pa- 
rameterization 



Tcoir) = Tco{ro){r/roy 



(A.12) 



while the mass density is calculated assuming that the 
disk is in hydrostatic equilibrium and vertically isothermal 
in the inner region at the midplane temperature T„j(r), 
through the relation 



p{r,z) = pair) 



(A.13) 



valid between the disk inner and outer radii i?i„ and Rout ■ 
The density on the disk midplane po{r) can by expressed 
as function of the surface mass density S(r) = So(r/ro)~^ 
through the relation 



Po{r) 



(A. 14) 



'27rh{r) 

Finally, the pressure scale h{r) is given by the relation 



hir) = 



/2r3kr,(r) 



(A.15) 



where Ti{r) is obtained solving the structure of a stellar- 
irradiated passive disk as described in DuUemond et 
al. (2001). 

Note that the emitting gas temperature Tco has been 
parametrized independently of the disk interior tempera- 
ture Ti, which governs the density structure of the disk. 
As pointed out by Dartois et al. (|2003p . the more optically 
thick CO transitions may by good tracer of the disk sur- 
face where the gas temperature is different from the disk 
interior. 

The resulting model, produces brightness maps for 
each frequency (see Fig |A.2p that can be compared with 
the observations presented in Fig. [2] and [3l 
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Fig. A. 2. ^•^CO J=l-0 simulated emission for a Keplerian disk inclined by 45°rotating around a central star of IMq. 
The disk outer radius is Rout — 100 AU with a surface density S(r) oc r~^'^. The left panel show the spatial integrated 
line profile. The right panel shows the spatially resolved maps at different velocity as labelled in figure. 
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